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The inverse energy cascade of two-dimensional (2D) turbulence is often represented phc- 
nomenologically by a Newtonian stress-strain relation with a 'negative eddy- viscosity'. 
Here we develop a fundamental approach to a turbulent constitutive law for the 2D in- 
verse cascade, based upon a convergent multi-scale gradient (MSG) expansion. To first 
order in gradients we find that the turbulent stress generated by small-scale eddies is 
proportional not to strain but instead to 'skew-strain,' i.e. the strain tensor rotated by 
45°. The skew-strain from a given scale of motion makes no contribution to energy flux 
across eddies at that scale, so that the inverse cascade cannot be strongly scale-local. 
We show that this conclusion extends a result of Kraichnan for spectral transfer and 
is due to absence of vortex-stretching in 2D. This 'weakly local' mechanism of inverse 
cascade requires a relative rotation between the principal directions of strain at differ- 
ent scales and we argue for this using both the dynamical equations of motion and also 
a heuristic model of 'thinning' of small-scale vortices by an imposed large-scale strain. 
Carrying out our expansion to second-order in gradients, we find two additional terms 
in the stress that can contribute to energy cascade. The first is a Newtonian stress with 
an 'eddy- viscosity' due to differential strain-rotation, and the second is a tensile stress 
exerted along vorticity contour-lines. The latter was anticipated by Kraichnan for a very 
special model situation of small-scale vortex wave-packets in a uniform strain field. We 
prove a proportionality in 2D between the mean rates of differential strain-rotation and 
of vorticity-gradient stretching, analogous to a similar relation of Betchov for 3D. Ac- 
cording to this result the second-order stresses will also contribute to inverse cascade 
when, as is plausible, vorticity contour-lines lengthen on average by turbulent advection. 



1. Introduction 

Almost forty years ago, Kraichnan (1967) predicted an inverse cascade of energy in two- 
dimensional (2D) incompressible fluid turbulence. This is perhaps one of the most intrigu- 
ing turbulent phenomena to occur in classical fluids. Kraichnan proposed an inertial range 
with a k~ 5 ' 3 power-law energy spectrum, just as in three dimensions (3D), but with a flux 
of energy from small-scales to large-scales rather than the reverse. Kraichnan's detailed 
predictions for steady-state forced 2D turbulence have been confirmed with increasing 
precision in a series of numerical simulations [Lilly (1971), Lilly (1972), Fyfe, Montgomery 
& Joyce (1977), Siggia & Aref (1981), Hossain, Matthaeus & Montgomery (1983), Frisch 
& Sulem (1984), Herring & McWilliams (1985), Maltrud & Vallis (1991), Boffetta, Celani 
& Vergassola (2000)] and laboratory experiments [Sommeria (1986), Paret & Tabeling 
(1998), Rutgers (1998), Rivera (2000)]. In fact, it can be rigorously proved that an in- 
verse cascade with constant (negative) flux of energy must occur in a forced 2D fluid, 
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if damping at low-wavenumbers keeps the energy finite in the high Reynolds number 
limit (Eyink (1996a)). Kraichnan's seminal concept of an 'inverse cascade' has since been 
fruitfully extended to other physical situations, such as inverse cascade of magnetic he- 
licity in 3D magnetohydrodynamic turbulence (Frisch et al. (1975)), of wave action in 
weak turbulence (Zakharov & Zaslavskii (1982); see also Zakharov (1967)) and of passive 
scalars in compressible fluid turbulence (Chertkov, Kolokolov & Vergassola (1998)). 

Attempts have often been made to account for the 2D inverse energy cascade phe- 
nomenon by a negative eddy-viscosity, either within analytical closure theories (Kraich- 
nan (1971a), Kraichnan (1971b), Kraichnan (1976)) or more phenomenologically (Starr 
(1968)). Such a description postulates a constitutive law for the turbulent stress pro- 
portional to the strain, = —IvtS^, with a viscosity coefficient vt < 0. However, an 
exact elimination of turbulent small-scales gives rise to a stress formula which is quite 
different: nonlocal in space, history-dependent and stochastic (Lindcnberg, West & Kot- 
talam (1987), Eyink (1996b)). Thus, any local and deterministic parameterization of the 
stress, such as by an eddy-viscosity, can be only an approximate representation at best. 
Nevertheless, such simplified constitutive relations can be quite useful to illuminate some 
of the basic physics of turbulent cascades and they are also important, of course, for use 
in practical numerical modelling schemes. 

In a previous paper (Eyink, submitted), hereafter referred to as (I), we developed a 
general approximation scheme for the turbulent stress, based upon a multi-scale gradient 
(MSG) expansion. We employed there the filtering approach to space-scale resolution 
in turbulence (Germano (1992)), which is also used in Large-Eddy Simulation (LES) 
modeling schemes (Meneveau & Katz (2000)). Within that framework we developed an 
expansion of the stress, first in contributions from different scales of motion and then 
in terms of space-gradients of the filtered velocity field. As a concrete application of the 
general scheme we considered in (I) the forward cascade of energy and helicity in 3D. 
In this paper we apply the same formalism to the 2D inverse energy cascade. In certain 
respects, the 2D theory is more difficult than 3D, because of certain peculiarities of the 
inverse cascade. We find that contributions to the stress from velocity-increments at 
sub-filter scales are much more important in 2D that in 3D. Also, terms second-order in 
space-gradients play a significant role in the 2D inverse cascade, whereas in 3D the terms 
first-order in gradients appear to suffice. Recognizing these facts has proved crucial to 
unravelling the physics of the 2D inverse energy cascade. 

However, 2D is simpler than 3D in respect of geometry. As we discussed in (I), the 
local energy flux is given in general by a scalar product 

77 = -S:r (1.1) 

where 5 is the filtered strain tensor and r is the deviatoric (i.e. traceless part of the) 
stress tensor t. The quantity 77 defined in (1.1) represents the rate of work done by the 
large-scale strain against the stress induced by the small-scales. In 3D, this expression 
involves three eigenvalues for each tensor, and also three Euler angles which specify the 
relative orientations of the tensor eigenframes. However, in 2D one has simply 

77 = -a(5T)cos(26), (1.2) 

where ±a are the two eigenvalues of 5, ±<5t/2 are the two eigenvalues of r, and 9 is 
the angle between the eigenframes of these tensors. We have taken 0^8^ ir/2 and 
(7, St ^ 0. Thus, the essence of the inverse energy cascade lies exactly in the tendency 
that ^ 9 < 7r/4. If arc the two eigenvectors of the deviatoric stress corresponding 
to the eigenvalues ±<5t/2, then there is a net tensile or expansive stress St/2 along the 
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ej|!" direction and a net contractile or compressive stress —5t/2 along the e_ direction. 
Therefore, when ^ 9 < 7r/4 holds, the stretching direction e+ of the strain is aligned 
primarily along the direction of net tensile stress, whereas the squeezing direction 
of the strain is aligned mainly along the direction of contractile stress. In that case, the 
stress cooperates with the strain rather than resists it, and negative work is done by the 
large-scales against the small-scales. 

Our primary objective in this work is to gain some understanding how this character- 
istic alignment comes about in 2D. In a negative- viscosity model, the stress is directly 
proportional to the strain or, equivalently, the alignment angle 6 = 0. This configuration 
leads to maximal inverse cascade but it is unlikely to occur uniformly throughout the 
flow. In fact, a main result of our work is that, to first-order in gradients and considering 
only the contribution to stress from scales of motion near the filter scale, the alignment 
is instead 9 — 7r/4 everywhere [Section 2.1.1]. We call such a stress law 'skew-Newtonian' 
and, from (1.2), it gives zero energy flux. Thus, to first-order in gradients, no energy flux 
can arise in 2D from strongly scale-local interactions, in agreement with a conclusion 
of Kraichnan (1971b). On the other hand, skew-Newtonian stress from smaller subscale 
modes can give rise to non- vanishing flux, since the stress is oriented at angle 7r/4 with 
respect to the strain at the same scale, not the large-scale strain 5 [Section 2.1.2]. We 
argue that the flux from skew-Newtonian stress produced by more distant subfilter scales 
is negative, on average, because of a relative rotation of the principal directions of strain 
at distinct scales. A plausible explanation for this characteristic rotation is advanced 
based on the exact equation for the rotation angle [Appendix A] and a heuristic model of 
'vortex-thinning' [Section 2.1.3]. Furthermore, two additional main mechanisms of inverse 
cascade are predicted by carrying our expansion to second-order in gradients: a Newto- 
nian stress with eddy-viscosity due to differential strain-rotation and a tensile stress 
directed along vorticity contour-lines [Section 2.2.1]. The latter effect was anticipated by 
Kraichnan (1976) [Appendix B] and it produces inverse cascade when vorticity-gradients 
are stretched by the large-scale strain. We derive an identity [Appendix C] that shows 
that, under the same condition, the eddy- viscosity due to differential strain-rotation is 
negative on average and produces inverse cascade. These mechanisms operate for stress 
produced by subfilter scales also, but more weakly the more distant in scale [Section 2.2.2]. 



2. The Multi-Scale Model in 2D 

In this section we shall develop for 2D the MSG expansion of the turbulent stress that 
was elaborated in general in (I). To keep our discussion as brief as possible, we shall 
refer to (I) for most of the technical details and only outline here the main points of the 
general scheme. We employ the standard 'filtering approach' (Germano (1992)), which is 
reviewed, for example, by Meneveau & Katz (2000). Thus, we filter the velocity field u 
with a kernel G at a selected length-scale I in order to define a 'large-scale' field u from 
scales > I and a complementary 'small-scale' field u' = u — u from scales < I. However, 
we further decompose the velocity field using test kernels r n (r) — £~ d r(r/£ n ) into 
contributions from length-scales > i n — \~ n l. The difference = — it^™ -1 ' 
then represents the velocity contribution from length-scales between £„_i and l n and 
yields a multi-scale decomposition 
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of the velocity field. In this paper we assume a scale-ratio A = 2. We also assume, for 
simplicity, that the kernels G and r are equal. Thus, the two filtered fields u and u = u^ 
at length I are equal and we need no longer keep the second as a distinct object. 

Since the filtered velocity fields u^ are smooth, they may be Taylor-expanded into a 
series of terms from mth-order gradients V" 1 ^™). Appropriate functionals of the velocity 
field may be expressed in this manner as a summation over both the gradient index m 
and the scale index n, which we call a multi-scale gradient (MSG) expansion. Among the 
most important quantities for which such a MSG representation may be developed is the 
turbulent stress tensor r. The latter quantity is defined mathematically as r = uu — uu. 
Physically, it gives the contribution of the small-scales to spatial transport of large-scale 
momentum and it is the quantity which requires 'closure' in the equation for the large- 
scale velocity u. It was proved in (I) that there is a convergent MSG expansion for the 
stress tensor, under realistic conditions for turbulent cascades. 

We should remark that two related but distinct approximations for the subscale stress 
were developed in (I). The first (I, Section 3) was a systematic expansion, which we 
shall refer to simply as the MSG expansion. This is a doubly-infinite series in orders of 
space-gradients and in scales of the velocity field, which converges to the exact subscale 
stress. However, as discussed in (I), the rate of convergence of the expansion in order m 
of space-gradients is apt to be quite slow as the scale-index n is increased. To obtain a 
more rapidly convergent gradient-expansion in the small-scales, we developed also a more 
approximate method (I, Section 4). In this modified approach the small-scale stress was 
estimated from velocity-increments for separation vectors in a certain subset for which 
the gradient-expansion is rapidly convergent, at all scales. The hypothesis underlying 
this approximation is that the stress due to velocity-increments for separation vectors 
from all subregions is similar and can be estimated, to a good approximation, by the 
stress arising from the distinguished subset. We referred to this modified expansion in 
(I) as the Coherent- Subregions Approximation (CSA), or the CSA-MSG expansion. It is 
guaranteed to converge rapidly, but its accuracy depends upon the quality of the basic 
hypothesis. The latter seems plausible but should be subjected to empirical tests. 

As we shall see below, it is more important to consider the contributions of subfilter 
scales in the 2D inverse energy cascade than it is in the 3D forward cascade. Therefore, the 
rapid convergence of the CSA-MSG expansion at small-scales makes it more practical 
than the systematic expansion for 2D, and only the former will be considered here. 
However, given the close formal relation between them, most of our qualitative, physical 
discussion below can be carried over, with some minor changes, to the systematic MSG 
expansion, and it is only for the purpose of quantitative comparisons that the CSA 
expansion is to be preferred. To describe this approximation scheme it is necessary to 
decompose the turbulent stress as t = g — u'u', where we refer to g as the 'systematic' 
contribution to the stress and to —u'u' as the 'fluctuation' contribution. For further 
discussion of these two terms and for mathematical formulas, see (I; 2.13-14). In terms 
of these two quantities, the general CSA-MSG expression for the stress in any dimension 
d was given in (I), to nth-order in scale index and mth-order in gradients, as: 

n n 

rl"' ro) = Q [ * Um) - u'} kUm) u'} k ' Um) . (2.2) 

fe=0 fe,fe'=0 

Using the results for m — 2 as illustration, as in (I), we have 

ra,(2) cf „ 2 duW duw cf „ 4 d* u w 

Q* = —l-'k^. q \- TTJTj — oT * 



dxi dxi 2d(d + 2) dxidx m dxfix 
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Ad{d + 2) 
and 

M ;W,(2) = _1 c^itAuM (2.4) 

m 

The coefficients C p in this model for p = 2,4, ... represent the partial pth-moments of 
the filter- kernel G over a spherical shell of separation vectors of length ss ^ , corrected 
by a multiplicative factor of iVfe = 2 kd to compensate for the decreasing volume of those 
shells with increasing k. Explicit expressions were given for these coefficients with a 
Gaussian filter, in (I), Appendix C f. Notice that, with the volume-corrected coefficients 
used here, the 'fluctuation' terms in (2.2) are decreased relative to the 'coherent' terms 
by the factors l/^/NuNk' . These were proposed in (I) as a consequence of a central limit 
theorem argument for the averages over volume that define the 'fluctuation' velocities in 
(2.4). Because of this, those terms are expected for larger k to be negligible relative to 
the 'systematic' contributions in (2.2). 

This brief synopsis provides enough background on the MSG expansion for our appli- 
cation in this paper to the 2D inverse cascade. For mathematical derivations and more 
extensive physical discussion, see (I). 

2.1. The First- Order Model 

To begin our discussion of the 2D energy cascade, we shall consider the CSA-MSG 
expansion of the stress developed to first-order in velocity-gradients. According to the 
general formula in equations (2.2)-(2.4), the expansion of the stress then contains only 
the 'coherent' part q, since the 'fluctuation' velocity u' vanishes to first-order. Thus, in 
any space dimension d, the expansion is given to this order by 



with 



consisting of just the first term in (2.3). See also (I; 5.2). Terms for large values of k 
become negligibly small (UV scale- locality) , so that the limit as n — > oo exists. For a 
monofractal velocity field with Holder exponent everywhere 1/3 — as expected in the 2D 
inverse cascade (Paret & Tabeling (1998), Yakhot (1999), Boffetta, Celani & Vergassola 
(2000))— the fcth term in (2.6) scales as - f k /3 (Eyink (2005)and I). 

We now specialize the model to 2D, using the standard formula for a velocity-gradient 
(deformation) matrix in 2D, 

dui 1 

which relates it the symmetric, traceless strain matrix S and (pseudo)scalar vorticity u>. 

f The expressions involve incomplete Gamma functions. For the case d = 2 relevant here, 

these become, for p = 2m, 7 x) = 7(1 + m,x) = ml [l - (l + x + ^ H + e^j , 

in terms of elementary functions. See Abramowitz & Stegun (1964), formulas 6.5.2 and 6.5.13 



6 



Gregory L. Eyink 



Here is the antisymmetric Levi-Civita tensor in 2D. Substituting (2.7) into (2.6) yields 

^^^{W+^ + ^M* 1 ! 8 } (2-8) 

where we have defined the skew-strain matrix as = Sifcefcjf . In terms of matrix arrays 

Thus, the skew-strain is also symmetric and traceless. It is easy to see that the strain 
and skew-strain are orthogonal in the standard matrix inner-product S: S = (and 
hence the prefix 'skew'). The various terms that appear in (2.8) are the same as those 
in equation (I; 5.4) for 3D and have the same physical interpretations. Note, however, a 

\k] \k] 

principal difference with 3D is the absence of terms proportional to uj\ or . Since the 
only component of vorticity is perpendicular to the plane of motion, no stress can be 
directed along vortex- lines in 2D. 

2.1.1. The Strong UV-Local Terms 

It is interesting to consider separately the first term in (2.5), for fc = 0, since it 
corresponds to the stress contribution from filter-scale velocity-increments. Thus, we refer 
to this as the strongly UV-local contribution. It is the only summand in the formula (2.5) 
which is closed in terms of the filtered velocity u = u^°\ In fact, this term corresponds 
just to the well-known Nonlinear Model for the turbulent stress (Meneveau & Katz 
(2000)), as discussed at length in (I). 

The most important observation about the strongly UV-local term in 2D is that it 
gives zero energy flux, pointwise in space. This is obvious for the term proportional to 
|cJ| 2 , since it is a pressure contribution. Furthermore, the first term is proportional to 
S 2 = a 2 1 in 2D, where / is the identity matrix, and is thus also a pressure contribution. 
Here we have used the Cayley-Hamilton theorem and the fact that the strain matrix in 
2D has two eigenvalues ±a of equal magnitude but opposite sign. Therefore, the first term 
contributes also zero flux. The term in (2.8) proportional to the skew-strain is deviatoric 
but it does not contribute to energy flux, by the orthogonality mentioned earlier. We can 
thus conclude that there is no energy flux anywhere in space arising from the strongly 
UV-local interactions, to first-order in velocity-gradients. 

This conclusion agrees with a result of Kraichnan (1971b), who showed that energy 
cascade in 2D cannot be strongly scale-local. It is worthwhile to summarize his demon- 
stration, which is based on the detailed conservation of energy and enstrophy in Fourier 
space. Let T(k,p,q) represent the energy transfer into wavenumbcr magnitude k from 
all triads of wavenumbers with magnitudes k,p, q. A measure of the scale-locality of the 
triad is provided by the parameter 

v = log 2 (fc mod /fc min ) ^ 0, 

where fc m in, fcmcd, W are the minimum, median, and maximum wavenumber magni- 
tudes, respectively, from the triad k,p, q. Intuitively, this quantity represents the 'number 
of cascade steps' between the minimum and median wavenumber. Note that fc max < 2fc me d 
by wavenumber addition, so that log 2 (fc max //c m i n ) ^ v+ 1. Thus, the parameter v unam- 
bigously measures the ratio of scales involved in the triadic interaction. In these terms, 

t This differs slightly from the general definition given in (I), which would lead us in 2D to 

term as 'skew-strain' instead the product lo 1 '5 . This slight difference in terminology should 
cause no difficulty. 
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nonlocal:): interactions correspond to those triads with v 3> 1 and strongly scale-local 
ones to those with 1/ « 1. Kraichnan (1971b) noted that in 2D the transfer function 
satisfies both 

T(k,p, q) + T{ Pl q, k) + T(q, k,p) = (2.10) 
as a consequence of energy conservation, and 

k 2 T(k,p, q) + p 2 T(p, q, k) + q 2 T(q, k,p) = (2.11) 

by conservation of cnstrophy. Multiplying through (2.10) by q 2 and subtracting from 
(2.11) gives 

(k 2 - q 2 )T(k,p, q) + (p 2 - q 2 )T{p, q, k) = 0. (2.12) 

Thus, if k 7^ p = q, then T(k,p,p) = 0, and substituting back into (2.10) gives also 
T(p,p,k) = T(p,k,p) = 0. Hence, there is zero transfer, if any two wavenumbers have 
equal magnitudes, and, in particular, if v = 0. However, it is very plausible to expect 
that the transfer function will be continuous in the wavenumber magnitude. In that 
case, transfer will be vanishingly small also in the limit that v« 1. Kraichnan (1971b) 
obtained more quantitative results using his analytical Test-Field-Model (TFM) closure. 
He found (see his Figure 2) that roughly 90% of the energy flux comes from triads with 
v > 1, 70% with v ^ 2, and 60% with v ^ 3. To obtain 90% of the total energy flux in 
the TFM closure required including all triads with v ^ 5. Thus, the 2D energy cascade 
was predicted by the Kraichnan to be scale-local (cf. also the exact analysis in Eyink 
(2005)) but only weakly so. 

There is a fundamental relationship between our argument and Kraichnan's. This is 
best understood by recalling the form of the energy flux in 3D from the strongly local, 
first-order terms, equation (I; 5.11) [and see also Borue & Orszag (1998)]: 

i7(°> J ) = ^C 2 £ 2 |-Tr(S 3 ) + ^-aJ T Su;| (2.13) 

Both of these terms vanish in 2D, the second because of absence of vortex-stretching. As 
discussed in (I), the first term can also be related to vortex-stretching, at least in a space- 
average sense, by a relation of Betchov (1956). Of course, the lack of vortex-stretching in 
2D is also what underlies the conservation of enstrophy, used in Kraichnan's argument. 
The argument that we have given confirms Kraichnan's conclusion and extends it to be 
also pointwise in space. 



2.1.2. The Weakly UV-Local Terms 

From the preceding discussion we can see that any energy flux that arises to first order 
in gradients must be due to subfilter modes, with k > 0. Since the contribution from 
modes with k 3> 1 is also small, the flux comes primarily from the weakly local terms 
with k > 1. This contribution for each k ^ 1 can arise solely from the skew-strain term in 
the stress (2.8), since, by the same reasoning as above, the other two terms are isotropic 
stresses or pressures. The flux from modes at scale k is thus 

nl^^Lcfii^^rs 1 ^ (2.14) 

\ Note that this definition makes no distinction between ultraviolet (UV) and infrared (IR) 
nonlocal interactions, as in Eyink (2005). 
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This can be rewritten in more intuitive fashion using 'polar coordinates' for strain ma- 
trices: 

cos(2a) sin(2a) \ ~q _ f — sin(2a) cos(2a) \ , , 

sin(2a) -cos(2a) j 1 b ~ ° \ cos(2a) sin(2a) J' ( > 

Here a = \S\/y/2 is the (positive) strain eigenvalue and a = (1/2) arctan(Si2/Sn) is the 
angle made by the frame of strain eigenvectors , with a fixed orthogonal frame. 
Note, incidentally, that the skew-strain is obtained by rotating the frame of the strain 
by 7r/4 radians. By choosing appropriately among the two unit eigenvectors ±e+\ one 
can always ensure that ^ \a\ < n/2. Thus, from (2.14), (2.15), 

ijW'W = cf ^°V fc U fc J sin[J>(aW - a(°))], (2.16) 

a remarkably simple and compact result. 

The total flux from all scales k = 0, 1, n to first order in gradients is thus 

n 

k=l 

In order to achieve an inverse energy cascade, it must hold that the terms in the sum are 
negative on average, at least for k > 1. The sign of (2.16) is determined completely by 
the factor uj^ sin[2(aM — ct^ -*)], which depends upon the relative angle — . If we 
choose that ^ |cJ fe J — \ < tt/2, then this factor will be negative if the strain-frame at 
scale k lags the strain-frame at scale (a^ < a^) in regions where u>^ > 0, and leads 
(cJ fc l > a^) in regions where < 0. Under these conditions, the small-scale stress 
will cooperate with the large-scale strain and the latter will do negative work. Note that 
this is quite different from a 'negative-viscosity' mechanism, with a Newtonian stress 
proportional to strain S' fc ' . Instead, the crucial deviatoric component of the stress is of 

the form 7M5 , where 7^ = C^' '^w^ /2 has dimensions of a diffusion constant and 
may be termed 'skew-viscosity' 

We see that a contribution to inverse energy cascade at scale k requires an anti- 
correlation in the signs of u>^ and — . A plausible dynamical mechanism for 
this can be suggested, based upon the exact equation for the strain orientation angle: 

(a^fD^a^ = l w ( fc >Q( fc ) + V-K^ + ■■■, (2.18) 
8 

with D[ k ^ = d t + u^-V the advective derivative at scale k, 

QW = Ap^ = i[^«] 2 - 2[a«] 2 (2.19) 

the pressure hessian at scale k, and 

K (k) = i(VpWxV)M^ (2.20) 

a space transport term due to pressure forces. The • • • terms in (2.18) represent contri- 
butions from the turbulent stress due to modes at length-scales < Ik- See Appendix A 
for the derivation. According to (2.18), D t ~ £- 2/3 in a 2D inverse energy cascade 
range, so that the rotation-rate increases with increasing k. Since the pressure contri- 
bution V--K" ( fe ) is spatially-nonlocal and averages to zero, it can be treated as random 
noise. We shall likewise disregard the effect of subgrid terms • • • . Thus, the expected 
correlation will be created in strain-dominated regions with < 0, since a { -> there 
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rotates against the locality vorticity uj^ and faster for larger k. Since the flux (2.16) is 
proportional also to the strain magnitudes most of the cascade should occur 

in the strain regions where this counter-rotation occurs. 

2.1.3. A Heuristic Model 

A simple model problem may help to illuminate the basic mechanism of inverse energy 
cascade due to skew-strain. We shall consider the effect of a large-scale uniform straining 
field 

*°> = ( C 3. ) • <-> 

on a collection of small-scale vortices, each initially circular with support radius l n . 
The ith vortex in the assembly will be assumed to have initially a vorticity distribution 
wj"^(|r — T"j|) radially symmetric about its center r». Let us assume also that the small- 
scale vortices have each a single sign of vorticity, but with the net circulation of the array 
equal to zero: J2i I dr wj"'(r) = 0. Kraichnan (1976) considered a very similar model 
problem of "vortex-blobs" in order to illustrate the mechanism of asymptotic negative 
viscosities in his Test-Field Model closure. In Appendix B we review Kraichnan's "blob- 
model" and compare it with the present one. Suffice it to say here that it was crucial in 
Kraichnan's calculation to take vortex wave-packets with a very rapid sinusoidal variation 
in the vorticity. On the contrary, we require no such variation and a particular case of our 
model is an array of vortex patches with constant vorticity levels, each initially circular. 

The effect of the straining field on this set of small-scale vortices will be to deform 
them into elliptical form, elongated in the x-direction and thinned in the y-direction. 
Kida (1981) found this behavior in his exact solution of 2D Euler for an elliptical vortex 
patch in a uniform shear flow, whenever the strain a and vorticity level u satisfy \a/ui\ ^ 
(3 - V5)/[2(2 + 2V5) 1/2 } = 0.15. More generally, the same phenomenon appears in 
a rapid distortion limit for the case of a strong strain 3> maxi ||wj™'||oo- We can 
then ignore the self-evolution of the vortices and also their mutual interactions. This 
permits us to focus on a single vortex centered at r = o with radial vorticity profile 
cj[™](r). The vorticity level set initially at radius r is distorted into an ellipse whose 
equation is x 2 /a 2 + y 2 /b 2 = 1 with semimajor axis a = rexp[a^°H] and semiminor axis 
b = rcxp[— a(°H] at time t. 

The immediate result is that the energy of the small-scale vortex patch is reduced, 
as a consequence of conservation of circulation. The area inside each elliptical vorticity 
contour is preserved, but the length of the perimeter is increased. In order to keep the 
circulation constant, the circumferential velocity must decrease. For example, in the case 
of a circular vortex patch of constant vorticity-level w'"' with initial radius r = £ n , the 
patch evolves into an elliptical shape with circulation • nab = Avfr^a, where iJ™'(i) 
is the x-component of the circumferential velocity at time t. The second expression for 
circulation holds in the limit when a^ h » 1 and a » b, so that the perimeter of the 
elliptical vortex is approximately 4a and is nearly parallel to the x-axis. In that case, 

u [n] (t) = (tt/4)wM& = (7r/4)t>]*„exp[-(7 (0) t]. (2.22) 

A similar argument can be made for points interior and exterior to the vortex, with the 
result that the velocity is everywhere reduced by a common factor of exp[— a^t]. Thus, 
the kinetic energy of the vortex is also decreased. (Of course, a single vortex of definite 
sign would have infinite energy in the unbounded plane, due to divergence at infinity. 
Such far-field divergence is absent when considering the array of vortices with zero net 
circulation.) 
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(b) 

Figure 1 . Mechanism of vortex-thinning, (a) A large-scale strain field with stretching direction 
along the :r-axis and shrinking direction along the y-axis, and a small-scale vortex of positive 
circulation, initially circular; (b) The vortex elongated along the x-axis and thinned along the 
j/-axis, and its strain basis, rotated by —45° with respect to the large-scale strain. 

The energy lost by the collection of small-scale vortices is transferred to the large-scales. 
To see this, observe that the large-scale straining, in addition to reducing the velocity 
amplitude of the small-scale vortices, also rectifies the velocity direction. The velocity 
vector of the elongated vortices points almost entirely in the x-direction and very little in 
the y-direction. Indeed, the vorticity level curve initially at radius r for the profile u)^(r) 
now becomes, to leading order, a pair of straight, parallel lines y — ±6 = ±r exp[— a^°H]. 
Thus, the vorticity field approximates to u^(y,i) = exp[cr^ ^]) when a^H ^> 1. 

This is just the vorticity associated to a long, narrow shear layer with weakened velocity 

U N ( yj t ) = _ exp[-a (0) £]sign (y) / J n] (r) dr (2.23) 

Jo 

directed entirely along the x-axis. If the tensor product u^u'™' were integrated over 
space at the initial time, it would produce only a diagonal stress contribution: 

Tij(t = 0) = f dr uf\r)uf{r) = % • tt f" dr r\uf ] (r)\ 2 , (2.24) 

J vortex J 

where u^\r) = (1/r) J Q r puj^(p) dp is the tangential velocity around the vortex center. 
(Here we have integrated only over the body of the vortex, neglecting the contribution 
of more distant regions). However, after "rectification" there is a net stress component 

Tu(i) = 2a f dyuW(y,t)uW(y,t) 




(2.25) 



with all other components much smaller. This resultant stress reinforces the large-scale 
strain field, so that J dr II(r,t) = —Sij Ty < 0, and negative work is done by the large- 
scales against the small scales. 
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This simple model of inverse energy cascade illustrates the pattern of relative orien- 
tation of strain frames at distinct scales, which was discussed earlier. In fact, within the 
long, narrow shear layer created by thinning of a vortex there is a velocity-gradient (or 
deformation) tensor of the form 

DN M)= (° -^fo'*)), (2.26) 
with (du^/dy)(y,t) = — uM^{y,t). The corresponding strain matrix is 



sW(vt)-( ° -<" N (i/,*)/2 \ (2 27) 



which has eigenvectors 



(2.28) 



for o;["](y,t) > and with e!T , e_ reversed for u>^ n \y,t) < 0. See Figure 1, which 
illustrates the case of a vortex patch of positive (counterclockwise) circulation. The small- 
scale strain basis shown there is rotated relative to the large-scale strain basis by — 7r/4 
radians. If the vortex patch had had negative (clockwise) circulation, then the rotation 
would have been by +7r/4 radians instead. 

This same model also clarifies the origin of stress proportional to skew-strain in our 
general scheme. The skew-strain in such an elongated vortex is 

~S [n] (vt)-( Jn]M/2 ° ^ (2 29) 

Let us introduce a convenient space-average over the vortex of the form 

(o;H) = (2/& 2 ) fdy f dy' J n \y',t) = (2/e 2 n ) t dr f dp J n \p). (2.30) 
Jo Jo Jo Jo 

By the Cauchy-Schwartz inequality, (l/l n ) J Q C dr {f^ J n \p) dp} 2 ^ [(4/2)(u w )] 2 , 
and, furthermore, these two quantities will generally have a ratio within some specified 
bounds. It follows that, when a^°H 3> 1, 

eM n] )(S [n] ) - ( Tl f _l /2 ) , (2.31) 

where we have set m = T\i/t\. Thus, a (deviatoric) stress proportional to skew-strain 
arises naturally from a narrow shear layer produced by vortex-thinning. 

It is not completely obvious why small-scale vortices in a two-dimensional inverse 
cascade range should be elongated and thinned by large-scale strain. After all, in such 
a range <t^ <~ £~ 2 / 3 <c (4i)~ 2 ^ 3 ~ w'"' for I > £ n . Thus, the large-scale strain is weak 
compared with the vorticity at smaller scales, exactly the opposite as is assumed in the 
rapid distortion limit above. The vorticity at length-scale l n could be expected to respond 
more strongly to the larger strains a^ n 1 » w'"' from length-scales l n i <C l n . However, 
the large-scale strain, although relatively weak, is coordinated over large distances and is 
temporally coherent, with a typical lifetime of tg ~ £ 2 ^ 3 . By contrast, the strain from the 
smaller scales is random and uncoordinated and, furthermore, evolves on a much shorter 
time-scale ti , ~ (t n ') 2 ^ 3 - Thus, the small-scale vorticity can adjust very rapidly to the 
persistent large-scale strain, whereas it does not have time to adjust to the many, even 
more rapidly fluctuating strains from the still smaller scales. 
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Clearly, our simple model calculation does not reflect all of the complexities of the 
two-dimensional inverse cascade range. However, it gives a simple physical picture for the 
origin of stress proportional to skew-strain, which, we believe, is essentially the correct 
one. If the initial profiles of the vorticity, utf 1 ^ (r) for the ith circular vortex, are not 
constant in the radial distance r from the center, then vortex-thinning produces also 
large vorticity-gradients parallel to the compressing direction of the strain field. This 
second-order effect will be discussed in detail in the following section. 

2.2. The Second-Order Model 

We have seen that, unlike in 3D, the MSG expansion rl"' m ' ) to lowest order in space- 
gradients, m = 1, can only explain energy cascade if subfilter scales n 1 are considered. 
However, another possible mechanism may be terms of higher order in space-gradients 
with m ^ 2. To investigate this possibility, we develop in this section the 2D MSG ex- 
pansion to second-order in velocity-gradients. One can specialize the formulas (2. 3), (2. 4) 
to 2D, replacing the velocity derivative with strain and vorticity using (2.7). The result 
is: 



Jfc],(2) _ !^[fc]„2 



1 Wfc] 



16 



elk] q[k] 

il,m jl,m 



(3,0 



1 ^iid.u, 



+ 32^ 



4k] 
[fej 



^«|Va,W| a 



and 



,'[fe],(2) 



1 T,[k] 



ik] 



(2.32) 



(2.33) 



In the last term of (2.32) and also in (2.33) we have defined di = tijdj, the skew-gradient, 
which satisfies V-V = 0. This is the same operator that appears in the stream-function 
representation of a velocity m, — d^. Indeed, to derive the last term in (2.32) and the 
term in (2.33) we used the stream function ip^ and the Poisson equation — A -0^1 = uj^ 
in order to write Auf = —duo^. 



2.2.1. The Strongly UV -Local Terms 

As for the first-order expansion, we begin by considering just the strongly UV-local 
terms with n — 0. These give altogether (note that = Ci 0) for n = 0) 



,(0,2) _ 



1 



= -c 



(0)^2 



1 



+ 16 64 1 



+l[Cf -2(Cf f]t A (d^idju). 



(2.34) 



Let us consider the physical meaning of the various terms that appear. 

We have already considered the terms in the initial line of (2.34) that arise from 
first-order velocity-gradients and have shown that they give no contribution to energy 
flux. The second line is remarkably similar in appearance to the first. In fact, it is not 
hard to see that the first term proportional to Su^ m Sji jTn is an isotropic (pressure) term, 
by exactly mimicking the argument we gave earlier for the SuSji-tevm, separately for 
each value of the index m that is summed over. Of course, the final term proportional 
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to Sij\Vuj\ 2 is also a pressure. This leaves only the middle deviatoric term as possibly 
contributing to energy flux. It is interesting that this second-order term, 

T (o),p] = ( 1 / 16 ) C r(<V(Vw.V)3 ) (2.35) 

gives rise exactly to an 'eddy- viscosity'. To see this, it is easiest to use the 'polar coordi- 
nates' (2.15) for the strain and skew-strain. Together with the chain rule, this gives 

( Vu7- V)S = -2( VZZJ- Va)S + ( VZU- VA)S (2.36) 

with A = In <7. Of course, the second term proportional to skew-strain does not contribute 
to energy flux. Thus, up to such conservative terms, we obtain 

T (0),p] = _(i/ 8 ) C |<V(VuJ-Va)S+ _ = _ 2ut s^ (2.37) 

with v T — C4°^ 4 (VZU-Va)/16. This is a stress of Newtonian form, with an eddy- viscosity 
due to differential-rotation of the strain. Indeed, the eddy-viscosity coefficient vt is just 
proportional to the rate of rotation of strain along the direction of maximum increase of 
vorticity. 

The final term of (2.34) arises from the combination of the last term in (2.32) for 
ft = and the product of two terms in (2.33) for k = k' = 0. These together give a 
stress exerted along the direction parallel to the skew-gradient VZU. Equivalcntly, this 
stress is directed normal to the vorticity-gradicnt VtU, or along the level-sets or contour- 
lines of the vorticity. There are two opposing contributions, a tensile stress proportional 
to cf^ from (2.32) and a contractile stress proportional to (C^) 2 from (2.33). Which 
dominates could depend upon the choice of the filter kernel G. However, the concrete 
calculations in (I), Appendix C show that C = [cf ] - 2(Cf } )] 2 /32 > for a Gaussian 
kernel. We have also checked this to be true for a few other cases, e.g. an exponential 
filter G{r) = e~^/(2ir). At least for these choices we see that there is a tensile stress 
of strength C£ 4 \'VuJ\ 2 exerted by the small-scales along vorticity contour-lines. As we 
discuss in Appendix B of the present paper, this effect was anticipated in a calculation of 
Kraichnan (1976) for a simple model problem of a 2D vorticity wave-packet in a uniform 
strain field. This tensile stress along vorticity contours should be contrasted with the 
contractile stress — l 2 \u}\ 2 /2 exerted along vortex-lines in 3D, discussed in (I). 

The strongly UV-local terms in the stress thus can give a non-vanishing contribution 
to energy flux, at second-order in gradients. Indeed, 

n (o)\2\ = _c£ 4 (VuJ) t S(Vuj) - C'e 4 S:(VuJ-V)S (2.38) 

with C = [Cf ] - 2(Cf } ) 2 ]/32 and C = cf } /16. Using e T Se = -5 and (2.36), this can 
also be written as 

7li 0)j21 = Ce 4 {VLj) T S(VZJ) + 4C"£ 4 W 2 (VZJ-Va). (2.39) 

These are the only UV-local contributions to the energy flux at second-order. 

It is important to determine the sign of these terms, on average, to see whether they 
contribute to inverse cascade or direct cascade. In this respect, note that the first term 
in (2.39) is proportional to the negative of the rate of vorticity-gradient stretching by 
the large-scale strain. That is, if one considers the equation for the large-scale vorticity 
gradient, then it has the form 

D t \VoJ\ 2 = -2( W) T S(VcJ) + • • • , (2.40) 

where D t = dt and • • • denotes neglected terms due to the turbulent stress. Thus, 
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we see that the first term in (2.39) is negative (inverse cascade) precisely when vorticity- 
gradients are magnified, a connection already noted by Kraichnan (1976). Equivalently, 
inverse cascade requires the stretching direction of the strain field to tend to be 
parallel to contour-lines of the large-scale vorticity. Since we have already seen that 
the small-scales induce a tensile stress along the contour lines, the stress and strain 
cooperate in this alignment and negative work is done by the large scales against the 
small-scales. Equation (2.40) renders the required alignment plausible, since components 
of the vorticity-gradient parallel to the squeezing direction will tend to grow, according to 
this equation. Note that this tendency might be moderated somewhat by the small-scale 
stress terms which we have neglected in (2.40); cf. Van dcr Bos et al. (2002). 

The second term in (2.39) will be negative precisely when VcJ-Va < 0. This means 
that the strain frame must counter-rotate against vorticity changes, i.e. rotate clockwise 
moving in the direction of increasing vorticity. We do not have a direct dynamical ex- 
planation for this tendency, analogous to the one we gave above for vorticity-gradient 
stretching. On the other hand, we have found that there is a simple kinematic relation 
between the rates of differential strain-rotation and vorticity-gradient stretching in 2D: 

((VcJ) T S(VcJ)) = -(5:(V57-V)f) (2.41) 

or, equivalently, 

((VuJ) T S(VoJ)) =4(W 2 {VoJ-Va)). (2.42) 

Equation (2.41) [or (2.42)] is an exact 2D analogue of the 3D relation of Betchov (1956), 
and, like it, depends just on homogeneity and incompressibility of the velocity field. For 
a proof of the '2D Betchov relation' (2.41), see Appendix C. An important immediate 
consequence is that differential strain counter-rotation and vorticity-gradient stretching 
must occur together, on average, while differential strain co-rotation is associated with 
mean shrinking of vorticity-gradicntsf . 

The net energy flux from both terms in (2.39) is always negative (inverse cascade) 
when there is mean stretching of vorticity-gradients. Because of the Bctchov-like relation 
(2.42) it follows that (77j 0) ' [21 ) = (C+C')l 4 r, where F is the common average in (2.42) 
and 

C + C= ^Cf + l[cf - (Cf ) 2 ] > 0. (2.43) 
To prove inequality (2.43), note that c|°^ ^ by its definition. Furthermore, 

Cf ] = ( dr \r\ 2 G(r) < J ( dr G(r) ■ [ dr \r\ 4 G(r) < Jcf^ (2.44) 

J\r\>l y J\r\>l J\r\>l 

by the Cauchy-Schwartz inequality and normalization of G. This gives (2.43). Thus, for 
any filter, the net flux is negative when T < 0. The 2D Betchov relation furthermore 
gives the ratio of contribution to inverse cascade of the two terms in (2.39), as C/C. 
For a Gaussian filter this ratio is C/C' = (1/2) — (9/13)e -1 / 2 = 0.08, so that approxi- 
mately 92.6% of the mean of (2.39) comes from differential strain-rotation and 7.4% from 
vorticity-gradient stretching. 

f Because it is purely kinematic, the '2D Betchov relation' holds just as well in the enstrophy 
cascade range. As discussed in Eyink (2001) and Chen et al. (2003), forward enstrophy-flux 
is also associated with mean stretching of filtered vorticity-gradients. Thus, differential strain 
counter- rot at ion must also occur, on average, in the enstrophy cascade. 
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2.2.2. The Weakly UV-Local Terms 

The terms of the MSG expansion that are second-order in gradients contribute to 
energy flux already from the strongly UV-local modes. However, there are additional 
contributions at second-order from all the other subscale modes. Here we shall discuss 
the physical interpretation and significance of those. 

In fact, the various terms that appear in the expressions for the 2D model stress, 
(2.32) and (2.33), can be readily understood. The first term in (2.32), which is first-order 
in gradients, has already been discussed. In the next group of three 2nd-order terms, 
the first and last are both pressure contributions and do not contribute to energy flux. 
However, the middle term is deviatoric and can give rise to flux. Using the analogue of 
(2.36), 

(vJ fe ].V)S lfel = -2(Vo,[ fe ].VcJ fc ])S [fel + (VJ fc ].VA[ fc l)S [fcl , (2.45) 

this term can be split into two. The first is a Newtonian stress — 2i/^'s' fc ' with an eddy- 
viscosity coefficient 

4 fcl = (l/16)cf ] £ 4 k (VJ k] -VaW) (2.46) 
arising from differential strain- rotation at a length-scale ty.. The other term is of the 
'skew-Newtonian' form ' with skew-viscosity coefficient 

= (l/16)cf ] 4(VJ fe ]-VA[ fc l) (2.47) 

arising from differential strain-magnification at the same length-scale £k- Note that we 
have defined the logarithm of the strain eigenvalue or magnitude as A^ = IncrM. Since 
the velocity field in the inverse cascade range is monofractal with Holder exponent 1/3 
(Paret & Tabeling (1998), Yakhot (1999), Boffetta, Celani & Vergassola (2000)), it is 
not hard to see that both v$ and 7^ are of order 0(^ 3 ), as expected, f. The last 

term in (2.32) represents a tensile stress of mag nitudc +cf ] £ 4 k \VJ^\ 2 /32 exerted along 
contour- lines of the vorticity uj^ at length-scale Ik. 

There remains the 'fluctuation' contribution to the stress from (2.33). This can be 
best understood by summing over scales, to give u'J n ' 2 ^ = V^"*™^ with a fluctuation 
stream-function 

-Ik] 



k=0 v fc 

Note that the factor 1/y/Nk reflects the cancellations that are expected to occur in 
the space-integral for the contributions from modes at length-scale Ik (I). We see then, 
finally, that — V^*™^ VV'*™^ represents a contractile stress along the streamlines of xpi n \ 
This term opposes and, to some degree, cancels against the tensile stress terms in (2.32) 
exerted along the contour- lines of for k = 1, n. 

If the model stress is substituted into formula (1.1) for the flux, then there results: 

n (n,2) = J2^c l 2 k] £tJ k \S^:S lk] ) + lc{ k] £i(VJ^-Va^)(S^:S^) 



k=0 



-lc lk] £UVio^-vx^)(s^--S lk] ) + ^cfei (vJ fc ]) T s(°)(v^]) 

f To show this, use the formulas 2Va = Sn V5i 2 -5i 2 V5n ^ VA = S xl VSi 1+ Si 2 V5 12 ^ 

1+^12 '^ll+ , - , 12 

and the general estimates from Eyink (2005) and (I) 
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-(v4 n) ) T ^°\Vipi n) ) (2.49) 

This is our final CSA expansion result for the energy flux in 2D. In addition to the first- 
order term that appeared in (2.17), there are now second-order contributions arising 
from differential strain-rotation, differential strain-magnification, and vorticity-gradient 
stretching. The final term in (2.49) is expected to be much smaller than the others, 
because of the cancellations in space-averaging discussed above and additional cancel- 
lations in the sum over scales in (2.48). We expect that the first four terms contribute 
to inverse cascade. For small k, should be correlated to some degree with S^°\ so 
that the differential strain-rotation and vorticity-gradient stretching terms ought to have 
negative mean- values, for similar reasons as the corresponding k = terms discussed 
earlier. Like the first-order 'skew-Newtonian' term, the differential strain-magnification 
term vanishes for k = and can therefore be expected to be relatively smaller than 
the differential strain-rotation term. It is interesting to note that the latter has its sign 
determined by the quantity Vw'^'Va^ cos[2(a[ fe l — a*- -*)], closely related to the signed 
quantity uj^ sin[2(cJ fc l — a^)] that appears in the first-order term. The final term in 
(2.49) is the only one that we expect to have a positive mean (from vorticity-gradient 
stretching), but we have already argued that that term will be considerably smaller in 
magnitude. 

Note that the flux term in (2.49) from scale k gives at most a fraction of order 2~ 2fe / 3 to 
the net energy flux. This agrees with rigorous locality estimates (Eyink (2005)). However, 
the actual contribution is likely to be much smaller, since the correlations which produce 
the inverse energy cascade must weaken for k 3> 1. If the small-scales are isotropic, then 
the mean stress r' fc l from length-scale £k will satisfy: 

(r W )= l (T r[r[ fe ]])^, forfc»l. (2.50) 

In that case, if the large-scale strain S^ ' and the stress contribution are asymptot- 
ically independent for k 3> 1, then their mean contribution to the energy flux vanishes, 
since the deviatoric part of the stress is zero on average. The existence of an energy 
cascade requires a statistical correlation between the large-scale strain and the small- 
scale stress contributions from various scales, which becomes progressively weaker for 
increasing k. 

3. Discussion 

The theoretical expression that we have developed here for the turbulent stress yields 
many concrete testable predictions — both qualitative and quantitative — for the 2D in- 
verse energy cascade. Foremost, we predict that strain-frames at small scales should 
lag/lead those at large-scales, when the small-scale vorticity is positive/negative. A spa- 
tial analogue of this effect is that the strain eigenframes are predicted on average to rotate 
clockwise in the direction of increasing vorticity (differential counter-rotation). Likewise, 
we predict that there will be a positive mean rate of stretching of vorticity-gradients. 
More quantitatively, our final CSA-MSG formulas (2.2), (2.32), (2.33) for the stress and 
(2.49) for the flux may be compared in detail with results obtained from experiment or 
simulation. If the model survives such tests, then it may be a good point of departure 
for building a practical LES modelling scheme of the 2D inverse energy cascade. 

In our presentation above we have alluded only briefly to the dynamical mechanisms 
that can produce the various correlations and alignments that are postulated, e.g. based 
on the evolution equations of strain orientation-angles (2.18) and of vorticity-gradients 
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(2.40). Many of the mechanisms expected to operate in 2D have very close analogues in 
3D. Notice that vortex-stretching in 3D is a near relative of the vortex-thinning mecha- 
nism in 2D, which we discussed in section 2.1.3. However, the result is opposite, because 
the stretching process in 3D "spins up" the vortices and increases the kinetic energy 
in the small scales. Vorticity contour-lines in 2D can also be expected to lengthen on 
the basis of the same plausible statistical arguments that have been applied to vor- 
tex lines or other material lines in 3D (Taylor (1938), Batchelor (1952), Cocke (1969)). 
This already argues rather strongly for the stretching of vorticity-gradients in 2D incom- 
pressible turbulence and, via the Betchov-like relation (2.42), for differential rotation of 
strain counter to vorticity. On the other hand, in 3D rather more detailed understand- 
ing is available through simple Lagrangian models of the evolution of velocity-gradients 
(Vieillefosse (1982), Vieillefosse (1984), Cantwell (1992), Chcrtkov, Pumir & Shraiman 
(1999)). These phenomenological models have provided plausible dynamical explanations 
of the key alignments that are observed in DNS (Ashurst et al. (1987)) and experiment 
(Tao, Katz & Meneveau (2002)). Some of the difficulties in developing such understand- 
ing of the inverse energy cascade can be appreciated by considering the exact equations 
in 2D for Lagrangian time-derivatives of the velocity-gradients: 



Here we have considered separately the evolution of the vorticity and strain. We have 
also neglected the contribution of turbulent stresses to the evolution of filtered gradients, 
which may be an important feedback interaction with small-scales (Van der Bos et al. 
(2002)). The equations (3.1) lack the local self-stretching terms which play the key role in 
the analogous 3D equations. In fact, the Lagrangian evolution in (3.1) is entirely trivial 
except for the pressure hessian in the equation for the strain and the latter must play an 
essential role in the production of strain orientation alignments. More sophistication in 
the modeling of pressure is therefore likely to be required than in the 3D case (Vieillefosse 
(1982), Vieillefosse (1984), Cantwell (1992), Chcrtkov et al. (1999)). Furthermore, wc have 
seen that in the 2D inverse cascade, both higher-order gradient and multi-scale effects are 
important. Thus, it remains a challenge to develop a detailed dynamical understanding 
of the 2D inverse energy cascade. 
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Appendix A. Dynamical Equation for the Strain Orientation 

It is easy to see from the 'polar' representation (2.15) of the strain 5 that 2a = 
arctan(Si2/Sn). Since also a 2 = S^ 2 + S^, the Lagrangian derivative may be written as 

2a 2 D t a = S U (AS 12 ) - S 12 (D t Su) (A 1) 



We can evaluate the time rate of change from the equation (3.1) for the filtered strain, 
which neglects the contribution from turbulent stress. Substituting into (A 1) we get 

d 2 p 1 f du dv" 
dy dx , 



„ 9 — du 
2a 2 D t a = — 
dx 



dxdy 
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d 2 ^ 
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du dv\ d 2 p 1 / du /\- (A2) 

dx dy ) dxdy 4 \ dy dx J P 



where we used incompressibility in the last line and also to derive the next identity: 



du d 2 p dv d 2 p 
dy dx 2 dx dy 2 



du dv\ d 2 p d ( dp du dpdu\ d ( dp dv dp dv 

dx dy J dxdy dx \ dx dy dy dx J dy \ dx dy dy dx 

(A3) 

If (A 3) is used in (A 2) to eliminate the mixed partial derivative of pressure, then one 
obtains 



2~ 2 ~D ~ ^ f dv du 

° tCt 4 \dx dy 
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dp du 
dx dy 



dp du\ 1 d 



dp dv dp dv 
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(A 4) 



This last equation is equivalent to (2. 18), (2. 19), (2. 20) in the text. 



Appendix B. Vortex-Thinning and Negative Eddy- Viscosity 

Some while ago, Kraichnan proposed a physical mechanism to explain the origin of 
negative eddy-viscosities in 2D (see Kraichnan (1976), Section 5.) For this purpose he 
employed a simplified model of small-scale vortex wavepackets in a uniform, large-scale 
straining field. His aim was to understand the asymptotic effect of the small-scales on 
much larger scales, and not to give an account of the inverse energy cascade by scale- 
local interactions. Nevertheless, his ideas turn out to have much in common with our 
theory of the local cascade interactions. The model proposed by us in section 2.1.3 to 
explain the stress proportional to skew-strain is just a slight modification of Kraichnan's. 
Furthermore, his mechanism of 'negative viscosity' is essentially identical with that we 
found in the last term of our model stress, equation (2.34), which corresponds to a tensile 
stress along vorticity-contour lines. Here we shall review the calculation of Kraichnan 
(1976), in order to make more clear its relation to the present theory. 

Kraichnan's model of the small-scales was a Gaussian wave-packet of vorticity — called 
a 'blob' — or an 'assembly of uncorrelated blobs' (Kraichnan (1976)). The stream function 
of each blob was taken to have the form 

ij){x) = k~ 2 f(x) cos(fca;2) 
f(x) = eM-l{A+A)/D 2 ) 

where / is a Gaussian envelope function with a standard deviation <~ D that is modulated 
by an oscillating cosine with wavevector k pointing in the vertical e 2 -direction. A basic 
assumption is that kD 3> 1, so that the wavenumbcr of the packet can be regarded as 
nearly sharp. Calculating the small-scale velocity field from u = — V^», it is not hard to 
show that the leading component of the velocity is 

ui ~ k~ x f{x) sin(foE 2 ) (B2) 
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and of the vorticity-gradicnt is 

(Vu) 2 ~ -A:/(aj) sin(for 2 ) (B 3) 

asymptotically for kD ^> 1. Cf. Eq.(5.4) in Kraichnan (1976). Thus, the dominant com- 
ponent of the total stress T = J r = J uu is 

7~n=fc- 2 y dxiy da; 2 exp(-(a;? +x^)/L> 2 )sin 2 (fca;2) - 7rL> 2 /2fc 2 (B 4) 

for fcD 3> 1. That is, the dominant stress is positive, or tensile, and exerted along the 
horizontal direction e.\. This is perpendicular to the direction of the vorticity-gradicnt 
e 2 , or along the direction of the vorticity-contours. Thus, Kraichnan's 'blob model' leads 
to a result in agreement with our general conclusion. 

As a model of the large-scales, Kraichnan took a uniform strain field 

S=a( C °f® 3in ^M (B5) 
y sm(20) — cos(20) J 

with eigenvalues ±a and eigenframe oriented at an angle <f> with respect to the fixed 
coordinate frame. The stream function corresponding to this large-scale field is just 
V(x) = 7}X T Sx. Actually, Kraichnan kept the strain fixed with frame axes along the 
coordinate directions and instead rotated the wavenumber of the small-scale blob, as 
k = k[e\sm(f) + e 2 cos^]; Eq.(5.14) in Kraichnan (1976). This is physically more nat- 
ural, if one thinks of the small-scales as isotropic and the large-scales as having fixed 
anisotropy. However, it is mathematically equivalent to rotate the strain and it relates 
more easily to our analysis in the text. 

Kraichnan (1976) worked out in detail the energy balance for his simple two-scale 
model of the velocity field. The initial energy in the small-scales is 

E=\J \u\ 2 ~ (1/2) Tu ~ ^D 2 /4fc 2 (B 6) 

The effect of the straining field on the small-scale wavevector is to change its magnitude 

by 

dk 2 /dt = -2k T Sk = 2acos(20)fc 2 . (B 7) 

Thus, Kraichnan concluded that, to leading order, 

{dE/dt) t=0 = -7ra£> 2 cos(20)/2fc 2 . (B8) 

Cf. Eq.(5.8) in Kraichnan (1976) for the case that k = 1 and <j> = 0. This reduction in 
energy of the small-scale blob is a consequence of the transfer of its enstrophy to higher 
wavenumber. 

Kraichnan showed further that the energy budget was maintained by a deposit into 
the 'interaction energy' J v-u between the large-scale and small-scale velocity fields. 
In his calculation he rewrote the interaction energy as J Vuj, in terms of the large- 
scale strcamfunction V and small-scale- vorticity u, and considered the nonlinear self- 
interaction of the latter. He found that the small-scale vorticity field set up a secondary 
flow of four equal-strength vortices with alternating signs of circulation which, for = 0, 
reinforced the large-scale strain. In his own words: 

'If a small-scale motion has the form of a compact blob of vorticity, or an assembly 
of uncorrelated blobs, a steady straining will eventually draw a typical blob out into an 
elongated shape, with corresponding thinning and increase of typical wavenumber. The 
typical result will be a decrease of the kinetic energy of the small-scale motion and a 
corresponding reinforcement of the straining field....' 
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In this way, the energy loss from the small-scales that is observed in (B 8) can traced to 
a transfer of equal size into the interaction energy between large-scales and small-scales. 

This transfer can be shown to be equivalent to the scale-to-scale energy flux that we 
defined in (1.1). Indeed, using the fact that the large-scale velocity v is stationary and 
its velocity-gradient V« is uniform, we find that 



This is the area-integral of the quantity that appears in (1.1). We can use this expres- 
sion to easily verify the energy balance result from Kraichnan (1976). Substituting the 
stress from (B4) and the strain from (B5), one gets (dE/dt) t =o — — iraD 2 cos(20)/2fc 2 , 
in agreement with (B8). Note that the flux is negative and the small-scales lose energy 
only if |0| < 7r/4, whereas the flux is positive for ir/4 < \4>\ < tt/2. If one assumes that 
the angle <ft is random with an isotropic distribution and fc = ke 2 is fixed, then the 
average flux is {{dE /dt) t= o) ang — 0. Kraichnan (1976) had already noted this result and 
established its consistency with the mean growth of small-scale wavenumber magnitude 
or, equivalently, the mean stretching of small-scale vorticity-gradients. As we discussed 
around our equation (2.50), a mean energy flux under isotropic conditions requires statis- 
tical correlations between disparate scales. In Kraichnan's case where he assumed a very 
wide separation between the two scales of motion, it was realistic to assume negligible 
correlations and thus zero net transfer. However, this is an unrealistic assumption in the 
context of a local energy cascade, where the stress and strain in (1.1) get most of their 
contributions from adjacent scales (Eyink (2005)) and are highly correlated. 

It is interesting that the mechanism that Kraichnan identified as acting between distant 
scales can also be identified with several of the mechanisms that we have found in our 
analysis of local cascade interactions. Note that in Kraichnan's vortex-blob model 



to leading order for Dk ^> 1. Thus, we get agreement of (B8) with the fourth term in 
our formula Tjj"' 2 " 1 for the energy flux, equation (2.49), by taking £k = 1/fc there. The 
second term in (2.49) corresponding to differential strain rotation is zero in the vortex- 
blob model because the orientations of the strain-fields (both large-scale and small-scale) 
are uniform in space. However, we can equally well understand the energy flux in the 
blob model based upon the first term in (2.49) [the same as (2.14)] that corresponds to 
relative rotation of strain at disparate scales. Indeed, in the blob-model, the vorticity is 




(B9) 



(Vu) T S(Vw) = -ak 2 f 2 {x) si-n 2 {kx 2 ) 008(20) 



(B10) 



using (B3) and (B5). Integrated over space, this yields 




(Bll) 



ui(x) ~ —f(x) cos(fca;2) 



(B12) 



and the small-scale strain of the blob is 




(B13) 



to leading order. Thus it is not hard to calculate that 




(B14) 
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and integrated over space this gives also 

J ujS:S = -naD 2 cos(20)/2 (B 15) 

to leading order for Dk >• 1. Multiplying (B 15) by £| = 1/fc 2 , we get also agreement of 
(B8) with formula (2.14). It is intriguing to note that, before averaging over space, the 
two contributions from (B 10) and (B 14) are exactly out of phase. It is another simple 
exercise to verify that the third term in (2.49), from differential strain-magnification, is 
also non-zero in the blob model and gives a contribution of the same sort. 

Thus, it is clear that most of the terms in our CSA-MSG formula (2.49) are repre- 
sented in Kraichnan's blob model, in particular, the flux from skew-strain, from differ- 
ential strain-magnification, and from vorticity-gradient stretching. All of these can be 
produced by a single mechanism of 'vortex-thinning'. Our somewhat simpler model of 
vortex patches in section 2.1.3 also illustrates these same flux terms, except in the case of 
constant-vorticity patches, for which only the flux from skew-strain survives. The increase 
in wavenumber that was considered by Kraichnan in his blob model and the asymptotics 
Dk ^> 1 play no essential role in the skew-strain mechanism. Indeed, note that (B 12)- 
(B 15) for the blob model all have non-vanishing values at k = 0, whereas (B 10)-(B 11) 
tend to zero as k — > 0. 



Appendix C. 2D Betchov Relation 

For any incompressible or solenoidal field u in 2D we can define a corresponding 'strain' 
and 'vorticity' u/") via 



diii _ _ („) 1 



a,, -V J "- (C1) 

Observe our notation for partial derivative with respect to Xj, indicated by subscript j 
preceded by a comma. Likewise, we write d 2 Ui/dxjdx k — Uij kl etc. Using these notations 
and definitions, the first step in the derivation of the 2D Betchov relation is the following 
identity: 

di[S^fvi,kWj t ki] - dkiS^fv^kWj^ii] = ^lfvi,kiWj,ki ~ S^v^kkWj.u 

*• v ' S v ' 

® © 

i c(u) c( u ) 

+ SijfakWjM ~ %;k v i,kWj,u 

® ® (C2) 

Here u,v,w are all incompressible fields. The identity (C2) follows straightforwardly 
from the product of rule of differentiation. 

The terms labelled ©, ©, and © are easily calculated by the substitutions Vij — 
_ {l/^)e i3 uj {v \vi^ k = S>?) k - (l/2)€ i: jd k uj {v \ and v tikk = -eudiu)^, and similar 
substitutions for the field w. The term © becomes 

S&Wi,« = \S% ] {^!A^ W) +ti^ iv) ) . (C3) 

© becomes 

-Slfv^kkW^i = S^d^d^K (C 4) 
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and © becomes 

The last term © requires as an additional step to use the identity 

5 Si - \tijdku {u) = mjk = Ui,kj = - JeaSjW 1 '' 1 (C 6) 

to replace S["\ by Sj^-. Then using the same substitutions as for the other three terms, 
® becomes 

-Sgk^H = S^Sl^M + Isg^M^M + Je« 8 iW Wfl jW (») w W (C 7) 

We are now able to sum the contributions from all four terms, ©, ©, ©, and ©. In 
order to simplify the result, it is helpful to define the quantity 

T$' v) = diuWdjU,™ + ~S { £ k d k ^\ (C 8) 

Then the sum of the four terms yields, after some elementary algebra, 

-2d k [S\fv iik w jtll ] + 2d l [S$ ) v i , k w jM ] (C 9) 

Observe that the first term on the righthand of (C 9) is antisymmetric under the inter- 
change u <-» w. Thus, if we symmetrize (C9) in u and w, we get 

5i; ) ^r ) +Perm. = div[...] (CIO) 

where the sum on the lefthand side is over all six permutations of u,v,w and div [• • •] 
on the righthand side indicates a total divergence. It therefore follows that 

(S^T^+perm^O (Oil) 

where (• • •) denotes either an average over a homogeneous ensemble or a space- average 
with boundary conditions that permit integrations by parts (e.g. periodic). We call the 
relation (C 11) the generalized Betchov identity in 2D. Setting u = v = w gives 

(% ) T\r ) )=° ( ci2 ) 

where and T^' u ^ are now constructed from the field u alone. Equation (C 12) is 
equivalent to the 2D Betchov relation (2.41) or (2.42) stated in the text. 



2D Inverse Energy Cascade 



23 



REFERENCES 

Abramowitz, M. and Stegun, I. A. 1964 Handbook of Mathematical Functions. National 
Bureau of Standards. 

Ashurst, W. T., Kerstein, A. R., Kerr, R. M., and Gibson, C. H. 1987 Alignment of 
vorticity and scalar gradient with strain rate in simulated Navier-Stokes turbulence. Phys. 
Fluids 30 2343—2353. 

Batchelor, G. K. 1952 The effect of turbulence on material lines and surfaces. Proc. Roy. Soc. 
Land. A 213 349-366. 

Betchov, R. 1956 An inequality concerning the production of vorticity in isotropic turbulence. 
J. Fluid Mech. 1 497-504 (1956). 

Boffetta, G., Celani, A. AND Vergassola, M. 2000 Inverse energy cascade in two- 
dimensional turbulence: Deviations from Gaussian behavior Phys. Rev. E 61 R29-R32. 

Borue, V. AND Orszag, S. A. 1998 Local energy flux and subgrid-scale statistics in three- 
dimensional turbulence. J. Fluid Mech. 366 1-31. 

Cantwell, B. 1992 Exact solution of a restricted Euler equation for the velocity gradient 
tensor. Phys. Fluids A 4 782-793. 

Chen, S., Ecke, R. E., Eyink, G. L., Wang, X. and Xiao, Z. 2003 Physical mechanism of 
the two-dimensional enstrophy cascade. Phys. Rev. Lett. 91 214501. 

Chertkov, M., Kolokolov, I., and Vergassola, M. Inverse versus direct cascades in tur- 
bulent advection. Phys. Rev. Lett. 80 512-515. 

Chertkov, M., Pumir, A., and Shraiman, B. I. 1999 Lagrangian tetrad dynamics and the 
phenomenology of turbulence. Phys. Fluids 11 2394-2410. 

Cocke, W. J. 1969 Turbulent hydrodynamic line stretching: consequences of isotropy. Phys. 
Fluids 12 2488-2492. 

Eyink, G. L. 1996 Exact results for stationary turbulence in 2D: consequences of vorticity 

conservation. Physica D 91 97-142. 
Eyink, G. L. 1996 Turbulence noise. J. Stat. Phys. 83 955-1019. 

Eyink, G. L. 2001 Dissipation in turbulent solutions of 2D Euler equations. Nonhnearity 14 
787-802. 

Eyink, G. L. 2005 Locality of turbulent cascades. Physica D, 207 91-116. 
Eyink, G. L. 2006 Multi-scale gradient expansion of the turbulent stress tensor. J. Fluid Mech. 
to appear. 

Frisch, U., Pouquet, A., Leorat, J., and Mazure, A. 1975 Possibility of an inverse cascade 
of magnetic helicity in magnetohydrodynamic turbulence. J. Fluid Mech. 68 769-778. 

Frisch, U. and Sulem, P. L. 1984 Numerical simulation of the inverse cascade in two- 
dimensional turbulence. Phys. Fluids 27 1921-1923. 

Fyfe, D., Montgomery, D. and Joyce, G. 1977 Dissipative, forced turbulence in 2- 
dimensional magnetohydrodynamics. J. Plasma Phys. 17 369-398. 

Germano, M. 1992 Turbulence: the filtering approach. J. Fluid Mech. 238 325-336. 

Herring, J. R. and McWilliams, J. C. 1985 Comparison of direct numerical-simulation of 
two-dimensional turbulence with 2-point closure - the effects of intermittency. J. Fluid 
Mech. 153 229-242. 

Hossain, M., Matthaeus, W. H., and Montgomery, D. 1983 Long-time states of inverse 
cascades in the presence of a maximal length scale. J. Plasma Phys. 30 479-493. 

KlDA, S. 1981 Motion of an elliptic vortex in a uniform shear flow. J. Phys. Soc. Jap. 50 
3517-3520. 

Kraichnan, R. H. 1967 Inertial ranges in two-dimensional turbulence. Phys. Fluids 10 1417- 
1423. 

Kraichnan, R. H. 1971 An almost-Markovian Galilean-invariant turbulence model. J. Fluid 
Mech. 47 513-524. 

Kraichnan, R. H. 1971 Inertial-range transfer in two- and three-dimensional turbulence. J. 

Fluid Mech. 47 525-535. 
Kraichnan, R. H. 1976 Eddy viscosity in two and three dimensions. J. Atmos. Sci. 33 1521- 

1536. 

Lilly, D. K. 1971 Numerical simulation of two-dimensional turbulence. Phys. Fluids Suppl. 12 
240-249. 



24 



Gregory L. Eyink 



Lilly, D. K. 1972 Numerical simulation studios of two-dimensional turbulence: I. models of 
statistically steady turbulence. Geophys. Fluid Dyn. 3 289-319. 

Lindenberg, K., West, B. J. and Kottalam, J. 1987 Fluctuations and dissipation in prob- 
lems of geophysical fluid dynamics, in: Irreversible Phenomena and Dynamical Systems 
Analysis in Geosciences, NATO ASI Ser. C, C. Nicolis and G. Nicolis, eds., D. Reidel 
Publ. Co. 

Maltrud, M E. and Vallis, G. K. 1991 Energy-spectra and coherent structures in forced 
2-dimensional and beta-plane turbulence. J. Fluid Mech. 228 321-342 (1991). 

Meneveau, C. and Katz, J. 2000 Scale-invariance and turbulence models for large-eddy sim- 
ulation. Ann. Rev. Fluid Mech. 32 1-32. 

Paret, J. AND Tabeling, P. 1998 Intermittency in the two-dimensional inverse cascade of 
energy: Experimental observations. Phys. Fluids 10 3126-3136. 

Rivera, M. K. 2000 The inverse energy cascade of two-dimensional turbulence. Ph. D. thesis, 
University of Pittsburgh. 

Rutgers, M. A. 1998 Forced 2D turbulence: experimental evidence of simultaneous inverse 
energy and forward enstrophy cascades. Phys. Rev. Lett. 81 2244-2247. 

SlGGlA, E. D. and Aref, H. 1981 Point-vortex simulation of the inverse energy cascade in 
2-dimension turbulence. Phys. Fluids 24 171-173. 

Sommeria, J. 1986 Experimental study of the two-dimensional inverse energy cascade in a 
square box. J. Fluid Mech. 170 139-168. 

Starr, V. P. 1968 Physics of Negative Viscosity Phenomena McGraw-Hill. 

Tao, B., Katz, J., and Meneveau, C. 2002 Statistical geometry of subgrid-scale stresses 
determined from holographic particle image velocimetry measurements. J. Fluid Mech. 
467 35-78. 

Taylor, G. I. 1938 Production and dissipation of vorticity in a turbulent fluid. Proc. Roy. Soc. 
Lond. A 164 15-23. 

Van der Bos, F., Tao, B., Meneveau, C, and Katz, J. 2002 Effects of small-scale tur- 
bulent motions on the filtered velocity gradient tensor as deduced from holographic PIV 
measurements. Phys. Fluids 14 2456-2474. 

Vieillefosse, P. 1982 Local interaction between vorticity and shear in a perfect incompressible 
fluid. J. Physique (Paris) 43 837-842. 

Vieillefosse, P. 1984 Internal motion of a small element of fluid in an inviscid flow. Physica 
A 125 150-162. 

Yakhot, V. 1999 Two-dimensional turbulence in the inverse cascade range. Phys. Rev. E 60 
5544-5551. 

Zakharov, V. E. 1967 On the theory of surface waves. Ph. D. thesis, Budker Institute for 

Nuclear Physics, Novosibirsk, USSR. 
Zakharov, V. E. and Zaslavskii, M. M. 1982 Kinetic equation and Kolmogorov's spectra 

in a weak turbulent theory of wind waves. Izv. Akad. S. USSR, Atmos. Ocean Phys. 18 

747-753. 



